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ABSTRACT 

We use the 6C** sample to investigate the co-moving space density of powerful, steep- 
spectrum radio sources. This sample, consisting of 68 objects, has virtually complete 
X-band photometry and spectroscopic redshifts for 32 per cent of the sources. In order 
to find its complete redshift distribution, we develop a method of redshift estimation 
based on the K — z diagram of the 3CRR, 6CE, 6C* and 7CRS radio galaxies. Based 
on this method, we derive redshift probability density functions for all the optically 
identified sources in the 6C** sample. Using a combination of spectroscopic and es- 
timated redshifts, we select the most radio luminous sources in the sample. Their 
redshift distribution is then compared with the predictions of the radio luminosity 
function of Jarvis et al. We find that, within the uncertainties associated with the 
estimation method, the data are consistent with a constant co-moving space density 
of steep-spectrum radio sources beyond z > 2.5, and rule out a steep decline. 

Key words: galaxies: active - galaxies: evolution - radio continuum: galaxies - 
quasars: general. 



1 INTRODUCTION 

Powerful radio sources, such as radio galaxies and quasars, 
trace the most massive galaxies (Jarvis et al. 2001a; De 
Breuck et al. 2002; Willott et al. 2003; Zirm, Dickinson & 
Dey 2003) and are associated with the most massive black 
holes (Dunlop et al. 2003; McLure et al. 2004; McLure & 
Jarvis 2004) in the Universe, at all cosmic epochs. Radio 
galaxies have been detected up to redshifts of just above 
five (van Breugel et al. 1999), and quasars up to and be- 
yond redshifts of six (Fan et al. 2003), leaving little time 
during which quasars and their host galaxies could form. 
This provides a challenge for hierarchical galaxy formation 
models, though some recent semi-analytic models are able 
to produce significant numbers of massive galaxies at high 
redshifts (e.g. Bower et al. 2006; Night et al. 2006). Con- 
straining the space-density of high-redshift radio sources is 
therefore important, as it has implications on the theories 
of galaxy and structure formation. Using radio sources is 
also advantageous in this respect because they are selected 
on the basis of their radio emission and are thus free of the 
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problems associated with optical selection methods, such as 
dust obscuration (cf. optically selected quasars). 

It is clear that the co-moving space densities of the 
rarest, most powerful radio sources were much higher around 
z ~ 2 than they are at present (Longair 1966), but the form 
of the evolution beyond that redshift is still a matter of 
debate. The question of a 'redshift cut-off' (Dunlop & Pea- 
cock 1990) in the radio source population has come under 
careful scrutiny in recent years (Jarvis & Rawlings 2000; 
Jarvis et al. 2001c). The current situation is one in which 
there is no compelling evidence of a high-redshift decline in 
low-frequency selected (i.e. predominantly steep-spectrum) 
radio sources, whereas there is evidence for a slight de- 
cline in radio-loud quasars from high-frequency selected (i.e. 
predominantly flat-spectrum) samples (Dunlop & Peacock 
1990, Shaver et al. 1996; Jarvis & Rawlings 2000; Wall et al. 
2005). 

In this study we focus on the low-frequency selected 
population (predominantly radio galaxies). To this effect, we 
use a new sample of radio sources drawn from the 151 MHz 
6C survey, which has been filtered with radio criteria de- 
signed to optimize the chances of finding radio galaxies at 
z > 4. This sample, namely 6C**, has been selected to be 
brighter than 0.5 Jy at 151 MHz. Additional selection crite- 
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ria have excluded all sources with radio spectral index be- 
tween 151 MHz and 1.4 GHz flatter than 1, or with radio 
angular size larger than 13 arcsec. The final sample consists 
of 68 objects over an area of sky of 0.421 sr, and is statisti- 
cally complete at an angular size limit of 9 < 11 arcsec. Full 
details of how the 6C** sample was selected can be found 
in Cruz et al. (2006, hereafter Paper I). 

The selection criteria just described are similar to those 
of the 6C* sample (Blundell et al. 1998; Jarvis et al. 
2001a, b), which was one of the samples used by Jarvis et 
al. (2001c) to constrain the co-moving space density of low- 
frequency selected radio sources. The 6C* sample was cru- 
cial in that study in sampling to high redshift (z ~ 4.4). 
The 6C** sample, being larger (cf. 0.13 sr) and deeper (cf. 
0.96 < 5*151 < 2.0 Jy) than 6C*, aims to improve on the 
small-number statistics limitation of this previous work, and 
ultimately to extend it to higher redshifts (z > 5) . 

Deep if- band imaging follow-up with UFTI/UIST on 
UKIRT, NIRI on Gemini and NIRC on Keck provided pho- 
tometry for all members of the 6C** sample (Paper I). Op- 
tical spectroscopy provided redshifts for 32 per cent of the 
sources (Paper I and references therein) . A summary of key 
observational information is given in Table 1. 

In this paper we describe a method of redshift estima- 
tion based on the K — z diagram of radio galaxies. This is 
presented in Section 2. In Section 3, we use the complete 
set of if-band magnitudes of the 6C** sample to estimate 
redshifts for all its optically identified members. These are 
compared to spectroscopic redshifts in Section 4, in order to 
assess the robustness of the method. The resulting estimated 
redshift distribution is discussed in Section 5. In Section 6 
we summarize the model radio luminosity function (RLF) of 
Jarvis et al. (2001c). This is the most relevant model to com- 
pare our data to, because it takes into account the selection 
effects of the 6C* sample. In Section 7 we compare the red- 
shift distribution (including spectroscopic redshifts) of the 
6C** sample with the model predictions, and discuss the 
evolution of the co-moving space density of the most radio 
luminous, low-frequency selected sources. Unless otherwise 
stated, we assume throughout that Ho = 70 km s _1 Mpc _1 , 
Qm = 0.3 and Qa = 0.7. The convention used for radio 
spectral index is S v oc v~ a , where S v is the flux-density at 
frequency v. 



2 REDSHIFT ESTIMATION METHOD 

Infrared-photometry provides a method of redshift estima- 
tion by utilising the tightness of the relation between K- 
band magnitude and redshift, which is characteristic of the 
near-infrared Hubble diagram of radio galaxies (Lilly & Lon- 
gair 1984; Eales et al. 1997; Jarvis et al. 2001a; De Breuck et 
al. 2002, Willott et al. 2003). The physical basis for the K- z 
relation is not well understood. At low redshifts, the Tf-band 
emission is dominated by the old stellar population in the 
host galaxy; at high redshifts, Jf-band samples rest-frame 
optical wavelengths, where the star formation history can 
have a significant effect. Non-stellar contamination to the 
if-band light, in the form of reddened quasar light and/or 
narrow emission lines, also contributes to the difficulty of in- 
terpreting the K — z diagram of radio galaxies, particularly 



at high redshifts (z > 3). Despite these caveats, the K — z 
diagram is still of interest as a tool for redshift estimation. 

Redshift estimates based on the K— z diagram have gen- 
erally been obtained by simple application of the empirical 
K — z relation (e.g. Dunlop & Peacock 1990). However, the 
significant amount of scatter around this relation requires 
the use of a more sophisticated method - one which takes 
into account all the available information in the diagram, 
and also which allows us to characterise the uncertainty on 
the output redshift estimates. With these requirements in 
mind the following approach is adopted: (i) we use Monte 
Carlo simulations to generate a statistical universe of syn- 
thetic realisations of the K — z diagram, based on a model 
of its underlying galaxy distribution, and (ii) we extract in- 
dividual photometric redshift probability density functions 
from this simulated population. 

2.1 The K - z Diagram for the 3CRR, 6CE, 6C* 
and 7CRS radio galaxies 

The most well defined K — z diagram for radio galaxies cur- 
rently available is the one obtained by Willott et al. (2003) 
from a combined dataset of the radio galaxies from the 
3CRR (Laing, Riley & Longair 1983), 6CE (Eales et al. 1997; 
Rawlings, Eales & Lacy 2001), 6C* (Jarvis et al. 2001a,b) 
and 7CRS (Lacy et al. 2000, Willott et al. 2003) flux-limited 
samples. It is based on a total of 204 radio galaxies with red- 
shifts ranging from 0.05 to 4.4, and its K — z relation is well 
fitted by a second-order polynomial between if-magnitudc 
and log 10 z (Willott et al. 2003): 

K(z) = 17.37 + 4.53 log 10 z - 0.31 (log 10 zf. (1) 

The main advantage of using this K— z diagram is that it has 
been obtained from completely identified samples with close 
to complete, or complete redshift information. This ensures 
the absence of significant biases in terms of sources with 
the weakest lines being missed because their redshifts are 
difficult to obtain. Another advantage is that these samples 
have been selected at a similar radio-frequency to 6C**, with 
progressively fainter flux-density limits. The brightest sam- 
ple is 3CRR selected at 178 MHz, with a flux-density limit of 
5*178 > 10.9 Jy (5isi > 12.4 Jy, assuming a spectral index of 
0.8); the faintest sample is 7CRS selected at 151 MHz, with 
a flux-density limit of 5isi > 0.5 Jy. The intermediate sam- 
ples are 6CE and 6C* selected at 151 MHz, with flux-density 
limits of 2.0 < S151 < 3.93 Jy and 0.96 < 5isi < 2.00 Jy, re- 
spectively. This results in a wide range in radio luminosity, 
which has made the investigation of the radio-luminosity de- 
pendence of the K — z relation possible in an unprecedented 
way. The correlation between if-band luminosity and radio 
luminosity has been one of the major worries with redshift 
estimates based on the K — z relation. 

Willott et al. (2003) found a statistically significant 
mean luminosity difference between the 3CRR and 7CRS 
radio galaxies of 0.55 mag in _R"-band, over all redshifts. How- 
ever, the 6C radio galaxies were found to differ on average 
from the 3C ones by only ~ 0.3 mag, which is much smaller 
than the value (~ 0.6 mag) reported previously (Eales et 
al. 1997). The mean luminosity difference between the 6C 
and 7C galaxies was not found to be significant. These re- 
sults are confirmed by McLure et al. (2004), who used HST 
data to study the host galaxy properties of a sample of radio 
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Table 1. Summary of the observational data on the 6C** sample. Column 1: Name of the 6C** source. Column 2: 151 MHz flux- 
density measurements in Jy from the 6C survey (Hales et al. 1988; Hales et al. 1990). Column 3: Radio spectral index evaluated 
between 151 MHz and 1.4 GHz from the 6C and NRAO VLA SKy Survey (NVSS; Condon et al. 1998) flux densities. Column 4: 
ii"-band magnitude within the angular aperture in arcsec given in brackets. Column 5: Rcdshift: '?' signifies that this value is uncertain, 
'nd' signifies that the source was observed but there were no emission lines and/or continuum detected. Column 6: Prominent emission 
line in the existing spectra, '?' signifies that the line identification is uncertain. Column 7: log 10 of the luminosity of the line listed 
in Column 6 (measured in units of W): '-' signifies that the data were inadequate to obtain a line luminosity, due to the absence of a 
spectrophotometric standard or non-photometric observing conditions. Column 8: Reference for the redshift of the source, Cea = Cruz 
et al. (2006), DBea = Dc Breuck et al. (2001), McC = McCarthy (1991), Vea = Vigotti et al. (1990), ASea = Allington-Smith et al. 
(1985), Bea = Brinkmann et al. (2000). 



galaxies at 0.4 < z < 0.6, which spans three decades in radio 
luminosity. They found mean luminosity differences in the 
i?-band of ~ 0.3 mag between the 3C and 6C samples, and 
~ 0.8 mag between the 3C and 7C samples 1 . These results 
are in good quantitative agreement with those of Willott et 
al. (2003). From both studies it follows that there is a corre- 
lation between radio luminosity and host luminosity within 
the 3CRR, 6CE and 7CRS samples. This has been inter- 
preted as suggestive of a correlation of both these proper- 
ties with black hole mass (Willott et al. 2003; McLure et al. 
2004). However, the weakness of these correlations means 
that the radio galaxies within the combined dataset used 
by Willott et al. (2003) follow essentially the same K — z 
relation. 



2.2 Parametric modelling of the K — z Diagram 

In this section we consider how to model the distribution of 
galaxies in the K — z diagram. We start from the combined 
dataset of radio galaxies from the 3CRR, 6CE, 6C* and 
7CRS samples 2 , and fit their distribution in log 10 z with a 
function of the form: 



n(z) — Aexp 



di (log 10 z) 



(2) 



where A and ai are free parameters. The histogram of the 
distribution in log 10 z along with the best-fitting function 
are shown in Fig. 1. We use a fifth order polynomial as 
the argument to the exponential, and find the best-fitting 
coefficients to be: 



A = 197.96; 



do = -0.39; 
ai = 1.17; 



a 2 = 1.00; a 4 = 1.47; 
a 3 = 1.83; a 5 = 0.38. 



For any given value of redshift z, we assume that n(z) of the 
sources follow a Gaussian distribution in K-band magnitude 
about a mean value k(z), given by the K — z relation (Eq. 1) 
at that redshift, i.e.: 



These magnitude differences were measured in the /-band and 
converted to the R-band assuming a constant R — I colour 
(McLure et al. 2004). 

2 These data can be obtained on-line at http: //www- 
astro. physics. ox. ac.uk/~cjw/kz/kz. html 




Figure 1. Histogram of the distribution in log 1() z of the ra- 
dio galaxies in the 3CRR, 6CE, 6C* and 7CRS samples. The 
bin width is A(log 10 z) = 0.25. The white regions correspond 
to 3CRR sources, the light shaded regions to 6CE sources, the 
intermediate shaded regions are the 7CRS sources, and the dark 
shaded regions the 6C* sources. The solid line represents the best- 
fitting function to the data - n(z) X A(log 10 z) — given by Eq. 2 
and the coefficients in the text. 



p(K\z) = 



n ( z ) 



exp 



[K-Hz)f 
2{a K ) 2 



(3) 



where K is the aperture- and emission-line corrected K- 
band magnitude measured for any given source (as described 
in Jarvis et al. 2001a; Willott et al. 2003) and gk , the dis- 
persion, is independent of redshift. This assumption is mo- 
tivated by the results of Willott et al. (2003), who evaluated 
the dispersion about the best-fitting K — z relation as a 
function of redshift. They found no significant increase up 
to 2 = 3, which is also in agreement with the results of Jarvis 
et al. (2001a). 

In order to constrain the free-parameter ok we fol- 
low the maximum likelihood formulation of Marshall et al. 
(1983). Defining S as —2 ln£, where C is the likelihood func- 
tion, the best-fit parameter is obtained by minimizing the 
value of S, which in our case is given by: 



The 6C** sample II: Redshift distribution 5 



S = -2Y,Hp{Ki\zi)] + 2 f f p(K\z)dK d(log 10 z). (4) 

i=i J 

In the first term the sum is over all N radio galaxies in 
the combined sample; the second term is the integral of the 
model distribution being tested and should give w 2N for 
good fits. The upper and lower limits of the integral are: 
10.0 < K < 21.0 and -1.3 < log 10 z < 1.0 (corresponding 
to 0.05 < z < 10). To find the best-fitting parameter we 
evaluated S over a range of 0.01 < Ok < 1- We found a 
minimum value of S for ok = 0.593 ± 0.02 mag. This value 
is in good quantitative agreement with the results of Willott 
et al. (2003), who found a scatter of 0.58 mag in the K — z 
relation at all redshifts up to z = 3. 



2.3 Redshift probability distributions from the 
simulated K — z diagram 

Adopting the model described in the previous section, 
such that p(K\z) dA"d(log 10 z) is the expected number of 
sources in the differential magnitude element dK and in 
the differential redshift element d(log 10 z), and using Pois- 
son probabilities, we use Monte Carlo simulations to gen- 
erate a large number of samples that mimic the combined 
3CRR/6CE/6C*/7CRS dataset. We combine all the sim- 
ulated datasets to construct a highly populated synthetic 
K — z diagram from which it is possible to extract photo- 
metric redshift probability density functions for any given 
value of K-magnitude. We extract these functions from a 
total of 10000 simulated samples, i.e. from a K — z diagram 
with ~ 2 million sources. 

The probability density functions p(z\K) can be ob- 
tained for any given source with 10 < K < 21 by taking 
the points along the horizontal band on the synthetic K — z 
diagram defined by [K - AK/2; K + AK/2], where K is the 
K-baxiA magnitude measured for that source and AK an 
appropriately small number, and by fitting the values of the 
relative frequency of each data point along this band. For 
any such source these values are best-fit by a log 10 -normal 
distribution of the random variable z with probability den- 
sity function: 



p(z\K) 



ln( 10) zV2tt a 2 



exp 



2<r 2 



(5) 



where /i and a are the mean and standard deviation of the 
distribution for the normal random variable log 10 z. The 
best-fitting estimate for z is thus defined as: 

zost = 10 M , (6) 

and the asymmetric 68% confidence interval as: 

10"-" < z cst < 10 M+CT . (7) 



the values for the best-fitting redshift estimates z est and 68% 
confidence intervals (C.I.), along with the p.d.f parameters a 
and n, and the if-band magnitudes used to extract them. In 
Fig. 2 we present the probability density functions p(z\K) for 
the 21 identified sources which have spectroscopic redshifts 4 
(Paper I and the references in Table 2). We note that there 
is in general good agreement between the redshift estimates 
and the spectroscopic redshifts. The notable exceptions are 
the quasars, which have systematically low estimated red- 
shifts because of the incorrect underlying assumption that 
the K-band light is dominated by starlight. 

Where possible we have used 7f-magnitudes measured 
in an 8-arcsec diameter aperture to extract the probability 
density functions. For sources which do not have if-band 
magnitudes measured in an 8-arcsec diameter aperture, due 
to the presence of a nearby object (Paper I), we have used 
those measured in 3- or 5-arcsec diameter apertures, and 
applied an empirical correction of —0.21 and —0.41 mag, 
respectively. These values have been derived from the me- 
dian difference between the small- (3 or 5 arcsec) and large- 
aperture (8 arcsec) magnitudes (presented in table 5 of Pa- 
per I). This choice of aperture and these corrections were 
designed to minimize the effects of the absence of an aper- 
ture correction in our analysis. 

We recall that the K—z diagram we use in our modelling 
is defined in terms of the aperture- and emission-line cor- 
rected Tf-band magnitudes. The aperture and emission-line 
corrections can only be obtained upon a priori knowledge of 
redshift and are thereby not available to us in this analy- 
sis. In the remainder of this section we will discuss how the 
absence of these corrections affects our redshift estimates. 



3.1 Aperture correction 

The aperture correction, e.g. as prescribed by Eales et al. 
(1997), involves converting the apparent angular size aper- 
ture magnitudes to standard 63.9 kpc metric apertures. This 
value was chosen because it corresponds to « 8 arcsec at 
z > 1 for a H = 50kms _1 Mpc -1 , fi M = 1 and Q A = 
cosmology. As it is immediately apparent from Fig. 3, aper- 
ture corrections are generally very small and for z > 0.6 
not strongly dependent on redshift. For the samples used 
in our modelling, at z > 0.6 magnitudes measured in 8- 
arcsec apertures have aperture corrections which are less 
than ±0.05 mag. For the magnitudes measured in 5- and 3- 
arcsec apertures, the values of the aperture correction are 
very similar to the empirical corrections we apply to our 
data. For these reasons, and for the high-redshift objects 
which we are most interested in, we consider the absence 
of an aperture correction on a source-to-source basis has a 
negligible effect on our redshift estimates. 



3 RESULTS 

Using K-band photometry from Paper I, we extract pho- 
tometric redshift probability density functions for all the 
identified sources in the 6C** sample 3 . In Table 2 we quote 

3 We exclude just two sources: 6C**0737+5618 and 
60**0935+4348. Both are not identified in K-band down to a 



3<r limiting magnitude of K ~ 21 in an 8-arcsec diameter aper- 
ture (Paper I). Therefore, we do not extract probability density 
functions for these sources, as the K — z diagram of radio galaxies 
is not well defined for K > 21 (see also the caption of Table 2). 
4 This excludes 6C**0935+4348, which has uncertain K-band 
identification and redshift. This source is discussed further in Pa- 
per I. 
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Redshift Redshift Redshift 

Figure 2. The redshift probability density functions p(z\K) for the 21 (identified) 6C** sources which have spectroscopic redshifts, 
normalised such that the area under the curve is unity. For each graph: the unshaded region corresponds to the 95% confidence in- 
terval; dotted lines show the location of the redshift estimate; dashed lines the location of the measured spectroscopic redshift. In 
some cases, where the spectroscopic redshift line is difficult to visualise, we add the label z Bpec . (a) 6C**0714+4616 (quasar); (b) 
6C**0726+4938; (c) 6C**0744+3702; (d) 6C**0746+5445; (e) 6C**0754+5019; (f) 6C**0810+4605; (g) 6C**0824+5344; 
(h) 6C**0832+5443; (i) 6C**0834+4129 (j) 6C**0854+3500; (k) 6C**0856+4313; (1) 6C**0903+4251; 
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Figure 2. continued: (m) 6C**0922+4933 (quasar); (n) 6C**0928+4203 (quasar); (o) 6C**0956+4735; (p) 6C**1009+4327; 
(q) 6C**1036+4721 (quasar); (r) 6C**1043+3714; (s) 6C**1045+4459; (t) 6C**1102+4329; (u) 6C**1132+3209. 



3.2 Emission line contribution 

Emission line contribution to the /("-band magnitudes can 
be significant at z > 2, when the strong Ha 6583, [O III] 
5007 and [O II] 3727 emission lines are redshifted into the 
if-band at 1.9 < z < 2.7, 2.9 < z < 3.7 and 4.2 < z < 5.2, re- 
spectively (Eales & Rawlings 1993, 1996; Jarvis et al. 2001a). 
The emission line correction, prescribed by Jarvis et al. 
(2001a), takes into account the [O II] emission- line-radio lu- 
minosity correlation of Willott (2001), and the emission-line 
flux ratios from McCarthy (1993) to determine the contri- 
bution to the if-band magnitude from all the emission lines 
at a given redshift. The emission line correction is therefore 
dependent both on radio luminosity and on redshift, but 
more strongly on redshift. This strong dependency is visi- 
ble in Fig. 4, where we plot the emission line correction as 



a function of redshift for the radio galaxies in the 3CRR, 
6CE, 6C* and 7CRS samples. However, it can be seen that, 
apart from one of the sources which has an emission line 
correction of 0.65 mag, all the sources in this dataset have 
emission line corrections which are smaller than the scat- 
ter in the K — z diagram. Therefore, we expect some of the 
most radio luminous sources in 6C**, which are at z > 2 to 
have an additional systematic uncertainty in their redshift 
estimates, due to emission line contamination to their K- 
band magnitudes. These will have redshift estimates, based 
on our method, which are biased towards lower values. This 
is a consequence of the contribution from bright emission 
lines, which make the objects brighter (in if -band). 
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Table 2. Redshift estimates for all the members of the 6C** sample. Column one lists the source names. Column two lists the i\-band 
magnitudes measured in, or corrected to an 8-arcscc diameter aperture (see Section 3). Columns three and four list the probability density 
function (Eq. 5) parameters fi and <r, respectively. Column five lists the redshift estimates as obtained through Eq. 6. Column six lists the 
68% confidence intervals as defined in Eq. 7. Column seven lists the spectroscopic redshifts (references as in Table 1). Columns eight and 
nine list the rest-frame 151 MHz radio luminosities (measured in units of W Hz -1 sr _1 and calculated assuming a power-law spectral 
index) based on the estimated and spectroscopic redshifts, respectively. Note: 6C**0737+5618 and 6C**0935+4348 are not identified 
in X-band down to a 3<r limiting magnitude of K ~ 21 mag in a 8-arcsec diameter aperture (Paper I). Therefore, we do not extract 
probability density functions for these sources. We give lower limits for their redshift estimates, although we caution that they have been 
determined by extrapolation of our method to K > 21 mag, where the K — z diagram is not well-defined. 
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Figure 3. The aperture corrections applied to the ii"-band magni- 
tudes measured in 8-arcsec (filled circles), 5-arcscc (filled boxes) 
and 3-arcsec (filled triangles) diameter apertures for the radio 
galaxies in the 3CRR, 6CE, 6C* and 7CRS samples as a function 
of redshift. The dashed lines correspond to the empirical correc- 
tions applied to the 6C** sources (see Section 3) of -0.21 mag 
(for magnitudes measured in 5-arcsec apertures) and -0.41 mag 
(for magnitudes measured in 3-arcsec apertures). 



4 COMPARISON WITH SPECTROSCOPIC 
REDSHIFTS 

In Fig. 5 we plot redshift estimates against spectroscopic 
redshifts for the 6C** sample. It is apparent that the ma- 
jority of sources lie under the line of equality between esti- 
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Figure 4. The emission line corrections applied to the if-band 
magnitudes (after aperture correction) of the radio galaxies in the 
3CRR (filled circles), 6CE (filled boxes), 6C* (stars) and 7CRS 
(filled triangles) samples as a function of redshift. 



mated and spectroscopic redshifts, implying a bias towards 
redshift under-estimation. This is not unexpected given the 
level of spectroscopic incompleteness of the sample. The sub- 
set of sources for which we have spectroscopic redshifts is 
biased towards sources which have strong emission lines and 
for which spectroscopy is easier to obtain. Depending on 
their redshift, these are more likely to be the sources in 
which emission line contamination to the if-band light is 
more significant. Indeed, if we exclude the quasars, it can be 
seen that the sources which show larger deviations from the 
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Spectroscopic Redshift 

Figure 5. Comparison of estimated and spectroscopic redshifts 
for 6C** sources. The error bars show the asymmetric 68% confi- 
dence interval about the best-fitting redshift estimate. The solid 
line is the z es t = Zspcc line. Triangles represent the sources with 
uncertain redshifts, boxes the sources with secure redshifts, and 
diamonds the quasars. The triangle with an upward arrow rep- 
resents the lower limit estimate for 6C**0935+4348 (see also the 
caption of Table 2) . 

equality line are in excess at z > 2. As we have discussed 
in Section 3.2, this is consistent with an unaccounted-for 
emission line contribution which systematically biases some 
of the sources towards lower redshifts. However, as it is 
clear from Fig. 2, the deviations are in all of the cases - 
6C**0746+5445 (z spcc = 2.156, Fig. 2d); 60**0832+5443 
(z spec = 3.341, Fig. 2h); 6C**0854+3500 (z Bpcc = 2.382, 
Fig. 2j) and6C**1045+4459 (z spec = 2.571, Fig. 2s) - within 
the 95% confidence interval upper-bound. 

Not surprisingly, the quasars as a whole show the largest 
deviations from the z CBt = z spcc line. The K — z diagram does 
not hold for these objects and therefore their K-bwaA mag- 
nitudes make poor redshift estimators. Quasars can have a 
strong contribution to their 7^-band light by the non-stellar 
continuum produced by the AGN. For this reason, any at- 
tempted estimates via the K — z diagram will be invariably 
much lower than the true redshift. This is particularly clear 
in Fig. 2 (a,m,q), where it can be seen that the estimates 
deviate by more than 2a, and in one case (60**0922+4216) 
more than 3<r, from the true redshift. We note, however, 
that the redshift of one of the quasars (60**0928+4203 
with z spoc = 1.664, Fig. 2n) is accurately estimated with 
our method. Possibly because its broad Mg II emission is 
the result of scattered light from the broad line region and 
the nucleus is heavily obscured (Paper I). 

Two sources have photometric redshifts which are 
significantly overestimated by our method. One is 
6C**1009+4327 (z spcc = 1.956, z ost = 3.089), which is a 
very faint source (K = 20.5 mag in a 3-arcsec aperture) 
with spectrum showing weak Lyman-a emission. The other 
is 60**0935+4348 (z spcc = 2.321, z cst > 4.0), which is one 
of the faintest sources in the sample, with K > 21.7 mag in 
a 3-arcsec aperture. This source is two magnitudes fainter 
than the mean K — z relation, which makes it a substantial 
K — z outlier (see Fig. 6). However, one caveat should be 
added here: the identification and redshift for this source are 




Redshift z 

Figure 6. K — z diagram of radio galaxies. The new spectro- 
scopically confirmed 6C** data are plotted along with the 3CRR, 
6CE, 6C* and 7CRS. The solid line is the K-z relation of Willott 
ct al. (2003). 6C**0935+4348, the major 6C** K - zoutlier, is 
represented by a box with an upward arraw, since this is also a 
ii"-magnitude lower limit. 

uncertain (Paper I). Sources like this are very unusual, but 
do exist. One of the 7CRS galaxies, 507.178 at z = 0.246, 
is also ~ 2 magnitudes fainter than the mean K — z relation 
(Willott et al. 2002b, 2003). 

For both 60**0935+4348 and 60**1009+4327, the dis- 
crepancy between if-band magnitude and redshift suggests 
that the stellar content is not consistent with that of a mas- 
sive host-galaxy. One possibility is that these sources are 
not yet fully formed (e.g. Jarvis, van Breukelen & Wilman 
2005). If, as a consequence of the 'youth-redshift' degeneracy 
described by Blundell & Rawlings (1999), we are preferen- 
tially observing young sources (<10 7 yr), then in some cases 
this could be the first instance of accretion activity in those 
sources. As such, this may trigger, or occur simultaneously 
with the initial star formation process in the host galaxy 
(e.g. Willott, Rawlings & Jarvis 2000; Willott et al. 2002a). 
If this is the case, it is perhaps not too surprising that their 
redshifts are overestimated by our method. As discussed by 
Willott et al. (2003) , the best-fitting K — z relation is close 
to the expected J^-magnitude evolution of a galaxy of local 
luminosity 3 L* which forms all of its stars in an instanta- 
neous burst at z/ = 10. 

However, 507.178, the 7CRS prominent outlier to the 
K — z relation, is at low redshift and therefore unlikely to 
be a forming galaxy. This suggests that the 'y° un g g a l ax y' 
hypothesis is not the only cause of underluminous outliers. 
The existence of a clear bright envelope to the K-z relation is 
presumably intimately related to the exponential cut-off in 
the galaxy Schechter function, whereas a tail to fainter mag- 
nitudes would be expected in any model in which powerful 
jets can, if only rarely, be associated with underluminous 
galaxies. 

Apart from the objects so far discussed, the major- 
ity of sources have been assigned with a reasonably well- 
constrained redshift estimate, in that the true redshift lies 
within the 68 per cent confidence interval about the best- 
fitting redshift estimate. Our method seems therefore to be 
fairly robust whenever emission- line and/or non-stellar con- 
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Figure 7. Estimated redshift distribution of the sources in the 
6C** sample (solid line). The bin width is Az = 0.5. The dark 
shadowed region is the distribution of sources in the 6C* sample 
(Jarvis et al. 2001b). The arrow above the 4.0 < z < 4.5 bin 
represents the lower limit redshift estimates for 6C**0737+5618 
and 6C**0935+4348. These are the two 6C** sources which are 
not detected in i\-band (see also the caption of Table 2). The light 
shadowed region shows the locus of the nine quasars of which we 
know about in the 6C** sample (listed in table 9 of Paper I). 



tributions to the K-b&nd light can be neglected. However, 
we note that this could be a selection effect resulting from 
the spectroscopic incompleteness of the sample. Among the 
subsample of objects for which we do not have spectro- 
scopic redshifts there could possibly be other sources like 
6C**1009+4327 and 6C**0935+4348, that have redshift es- 
timates significantly higher than their true redshifts. 

We now compare the near-infrared Hubble diagram of 
the 6C** radio galaxies (for which we have a spectroscopic 
redshift) with the 3CRR, 6CE, 6C* and 7CRS samples (see 
Fig. 6). We find that the 6C** objects follow the same 
K-z relation as the combined 3CRR/6CE/6C*/7CRS data, 
and have a rms dispersion of ~ 0.59 mag over all red- 
shifts if 60**0935+4348 is excluded. This is consistent with 
the dispersion value used in our modelling (Section 2.2). If 
60**0935+4348, the prominent K — z outlier, is included in 
this calculation, then the dispersion becomes ~ 0.75 mag. 



5 THE ESTIMATED REDSHIFT 
DISTRIBUTION 

The estimated redshift distribution for the entire 6C** sam- 
ple is presented in Fig 7. The median estimated redshift is 
z ~ 1.6. It is informative to compare these results with the 
redshift distribution of the 6C* sample (Jarvis et al. 2001b) 
which, as mentioned before, has been filtered in a very sim- 
ilar way. The 6C* sample has a median redshift of z ~ 1.9 
and a redshift distribution skewed towards z > 2. This is a 
direct result of the filtering criteria applied to it. 

Although the distribution of 60** sources does have 
a tail to high redshift, the great majority of sources 
are at z < 2, in clear contrast to what should be ex- 
pected. However, and in light of Section 4, we have to 
consider what is likely to be the major source of bias 



in our redshift estimation: the presence of quasars in 
the sample. These sources will be the ones which will 
be skewing the distribution the most towards lower red- 
shifts. The redshift estimates of the three quasars with un- 
resolved identifications (6C**0714+4616, 60**0922+4216 
and 6C**1036+4721), for which we have a spectro- 
scopic redshift, are under-estimated by an average fac- 
tor of 0.4. The major cause for concern is that there 
are quasars for which we do not have spectroscopic in- 
formation. On the basis of K-hand imaging alone, we 
assume that the further following sources are quasars 
(given that they have bright unresolved identifications; Pa- 
per I): 6C**0849+4658, 6C**1003+4827, 6C**1052+4349, 
6C**1056+5730 and 6C**1138+3803. This leaves us with 
nine quasars 5 , eight of which have redshift estimates in 
the range 0.5 < z < 1.0. These alone represent a significant 
source of the bias towards low-redshift which we witness in 
Fig 7. 

Removing the quasars from the distribution results in 
a median estimated redshift of z ~ 1.7. This is similar to 
what is found for 60* and much higher than the median red- 
shift of similar, unfiltered samples (e.g. the 7CRS sample, at 
the same flux-density limit, with z w 1.1). This result con- 
firms the efficiency of the filtering criteria, used for 60**, in 
excluding low-redshift sources. Furthermore, we recall that 
some sources may also have their redshifts systematically 
under-estimated, due to emission-line contamination to their 
K-band magnitudes (Sections 3.2 and 4). Thus, it is quite 
plausible that the real redshift distribution has a slightly 
higher median redshift. 

The fraction of quasars in the 60** sample (9 quasars 
out of 68 objects) is higher than that in the 60* sample. 
There are only two quasars out of the 29 60* objects (Jarvis 
et al. 2001b). Although the difference is barely significant 
given the small numbers involved, it is possibly influenced 
by small differences in the selection criteria, such as the 
tighter size constraint or, more importantly, the lower fre- 
quency spectral index cut. For 60*, the steep spectral index 
constrain is evaluated between 151 MHz and 4.85 GHz, thus 
excluding objects with prominent flat-spectrum cores, i.e. 
quasars. Because the 60** spectral index cut goes up to 
only 1.4 GHz, this effect is less pronounced. 



6 THE RLF MODEL OF STEEP-SPECTRUM 
RADIO SOURCES 

Jarvis et al. (2001c) investigated the radio luminosity func- 
tion for the most radio luminous low-frequency selected 
sources. The data used in their analysis were drawn from the 
complete 3CRR, 6CE and filtered 60* samples by imposing 
a lower limit in radio luminosity. Only the sources which lie 
in the top-decade in log 10 Li 5 i were considered 6 . Focusing 

5 These include 6C**0928+4203, which does not appear unre- 
solved in our near-infrared imaging, but shows broad emission 
lines in its spectrum (Paper I). We recall that this source has 
a redshift estimate which is in good agreement with its spectro- 
scopic redshift (Section 4). 

6 This meant: log 10 L151 > 27.63 in their cosmology I (Qm = 
1.0, tt A = 0.0, H = 50 km s" 1 Mpc" 1 ) and log 10 L 151 > 
27.90 in their cosmology II (fi M = 0.3, Q A = 0.7, H = 50 
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on only the most luminous sources provided the largest pos- 
sible baseline in redshift for the samples considered, while 
minimizing the role played (in the modelling procedure) by 
the intrinsic correlations between sample parameters, such 
as luminosity - spectral index and linear size - spectral index 
correlations, which are inherent to radio samples (Blundell, 
Rawlings & Willott 1999). 

The filtered 6C* sample represented the greatest ad- 
vance of this study over earlier work on the high-redshift 
RLF of low-frequency selected radio sources (Willott et al. 
2001). However, its use required that the effects of the selec- 
tion criteria, namely small angular size and steep radio spec- 
tral index, and in particular the fraction of sources which 
were missing from the survey, were taken into consideration 
in the modelling of the RLF. This led to a parameterisation 
which was separable in luminosity and redshift (as in Jarvis 
& Rawlings 2000 and in Willott et al. 2001), and also incor- 
porated distributions in radio spectral shape and linear size, 
i.e. 

p(L 151 ,z,ai,a 2 ,D) = p a x p L (L 15 i) x p(z) 

xp a (a 1 ,a 2 ) x pd(D), (8) 

where p is the normalising factor and a free parameter 
measured in units of Mpc~ 3 , pi(LiBi), p{z), p a (ai, a 2 ) and 
Pd(D) are dimensionless distribution functions in radio lu- 
minosity, redshift, spectral shape parameters and projected 
linear size, respectively. 

The full form of these distributions and their modelling 
are described in detail in Jarvis et al. (2001c). We now com- 
ment on the shape of the redshift distribution, for which 
Jarvis et al. (2001c) tested three different models. The model 
favoured by their maximum likelihood analysis (model C) 
uses a 1-tailed Gaussian to parameterise the low-redshift co- 
moving space density and a power-law distribution at high 
redshift, i.e 

where z is the 'break' redshift where the model switches 
from the low- to the high-redshift form, zi is the charac- 
teristic width of the half-Gaussian and rj is the power-law 
exponent describing the high-redshift co-moving space den- 
sity. The other models tested by Jarvis et al. (2001c) were: 
model A - parameterised as a single Gaussian distribution in 
redshift; and, model B - a 1-tailed Gaussian which becomes 
constant and equal to unity beyond the Gaussian peak, i.e. 
the same as model C with ij fixed at zero. The advantage 
of model C over the other models is that: (i) it breaks the 
symmetry between low- and high-redshift evolution which is 
forced by model A (and for which there is no physical basis) , 
and (ii) allows for freedom in the evolution at high-redshift 
which is not possible with model B. Model C, is therefore 
more useful in terms of assessing the form of the evolution 
of the co-moving space density at high redshift. 

Jarvis et al. (2001c) found a best-fitting power-law ex- 
ponent r\ = —0.06 (in their Cosmology II), implying a con- 
stant co-moving space density beyond a peak redshift of 

km s _1 Mpc -1 ), the two sets of cosmological parameters consid- 
ered in their analysis. 



z = 2.15 to an indeterminable redshift. A steep decline has 
been ruled out by their analysis at the ~ 4<r level, but the 
form of the evolution at high redshift remained unresolved 
with an uncertainty encompassing both moderate declines 
and continuing shallow inclines. 



7 COMPARISON WITH THE RLF OF 
STEEP-SPECTRUM RADIO SOURCES 

To investigate the co-moving space density of high-redshift 
steep-spectrum radio sources, we now compare the redshift 
estimates derived from the 6C** data with the RLF of low- 
frequency selected radio sources of Jarvis et al. (2001c). 

7.1 Selecting the most luminous 6C** sources 

To be able to compare the 6C** data with the RLF model of 
Jarvis et al. (2001c), we must select the most radio luminous 
sources in the sample in a fashion that is equivalent to that 
of Jarvis et al. (2001c). For that we use the definition of top- 
decade in cosmology II of Jarvis et al. (2001c) and translate 
its lower radio luminosity limit to log 10 L151 > 27.61 in 
the cosmology used in this paper (Q.m = 0.3, Q.a = 0.7, 
H = 70 km s^ 1 Mpc" 1 ). This means that, based strictly 
on the estimated L151 — z diagram (Fig 8: top-panel), 12 
of the 6C** sources are selected. However, if we consider 
also the spectroscopic redshifts available, i.e. by replacing 
estimated redshifts with spectroscopic ones where these are 
available (Fig 8: bottom-panel), then there are 17 sources in 
the top decade of radio luminosity. These sources are listed 
in Table 3, and include 6C**0737+5618. The inclusion of 
this source is discussed in more detail in Section 7.2.3. We 
note that two of the quasars with under-estimated redshifts 
(and luminosities) are now in the top-decade. A few of the 
quasars without spectra are also expected to move into the 
top-decade once their redshifts are known. 

7.2 Analysis and discussion 

We consider two ways of constructing the binned redshift 
distribution of 6C** sources in the top-decade of luminos- 
ity. In Section 7.2.2, we use spectroscopic redshifts when 
available, and otherwise the best-fitting redshift estimates, 
i.e. we consider the sources listed in Table 3 (see also Fig 8: 
bottom-panel). In Section 7.2.3, we take into account the 
redshift probability distribution of each source. Both red- 
shift distributions are then compared with the expected red- 
shift distribution given Jarvis et al. (2001c) model RLF. 

7.2.1 The predicted redshift distribution 

The predicted redshift distribution is obtained by integrat- 
ing the model C RLF of Jarvis et al. (2001c), normalized to 
the sky area and evaluated for the radio selection criteria of 
the 6C** sample, i.e. 

■■jjW = j j j j 'pC^isi, 2, ai, a 2 , D) 
xO(ii5i, z, 01, a 2 , D) 

x ^ d(log 10 L151) doi da 2 d(log 10 D) (10) 
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Source 



Estimated 



Spectroscopic 



CI. 



Estimated 



name 


Zest 


lo Sl0 L 151 


Zspec 


lo SlO L 151 




6C 0737+5618 


>4.0 


>28.18 








6C 0744+3702 


2.549 




2.992 


27.78 


G 


6C 0754+5019 


3.331 


27.97 


2.996 


27.85 


G 


6C 0801+4903 


2.722 


27.79 








6C 0824+5344 


2.507 




2.824 


27.71 


G 


6C**0832+5443 


2.240 




3.341 


27.69 


G 


6C**0902+3827 


2.246 


27.71 








6C**0909+4317 


1.881 


27.83 








60**0922+4216 


0.512 




1.750 


27.67 


Q 


6C**0925+4155 


3.060 


27.77 








6C**0928+4203 


1.736 


27.62 


1.664 




Q 


6C**0935+4348 


>4.0 


>28.38 


2.321 


27.75 


G 


6C**1009+4327 


3.089 


28.42 


1.956 


27.91 


G 


6C**1036+4721 


0.849 




1.758 


27.81 


Q 


6C**1045+4459 


1.569 




2.571 


27.66 


G 


6C**1050+5440 


2.773 


27.79 








6C**1102+4329 


2.732 


27.78 


2.734 


27.79 


G 


6C**1103+5352 


2.947 


28.21 








Table 3. Sources 


selected 


from the 6C 


** samp] 


c which corre- 



spond to our top-decade definition (see Section 7.1). The rest- 
frame 151 MHz luminosities are in units of W Hz -1 sr _1 and 
were calculated assuming a simple power-law spectral index. The 
third column corresponds to a top-decade based strictly on the es- 
timated redshift distribution. In this column the solid horizontal 
lines are used for sources which, based on their redshift estimates, 
fall below the lower luminosity limit. However, these sources 
are included in the top-decade if spectroscopic information (in 
columns 4-5) is taken into account. One source - 6C**0928+4203, 
is excluded. The classifications (G - radio galaxy; Q — quasar) 
are given only for sources which have optical spectroscopy (Paper 
I and references therein) and are based on these data. 



where p(L\^\, z, oi, a,2, D) is the complete radio luminos- 
ity function Jarvis et al. (2001c); fi(Li 5 i, z, oi, a,2, D) is 
the sky area available for the 6C** sample; and, (dV/dz) 
is the differential co-moving volume element. The factor 
f2(Li5i, z, oi, d2, D) assumes a value of 0.421 sr (the sky 
area of 6C**) or zero, depending on whether or not a source 
with a given a set of the parameters iisi, z, oi, a 2 and D 
meets the sample selecting criteria. The lower and upper 
limits of the integral are 27.61 < log 10 L151 < 28.61, —2.2 < 
01 < 1.0, -0.4 < a 2 < 0.2 and -0.3 < log 10 D < 4.0. The 
integral was evaluated numerically, between < z < 8, using 
a Monte Carlo method, over 10 7 random points uniformly 
distributed in the 4-dimensional parameter space. 

Because we use a different set of cosmological parame- 
ters than that of Jarvis et al. (2001c), we re-evaluated the lu- 
minosity function by using the relation from Peacock (1985): 



it \ dVl it \ dV 2 



(11) 



where pi and p2 are the luminosity functions for two differ- 
ent cosmologies, and, L\ and L2 are the luminosities derived 
from the flux density, redshift and proper motion distance 
in each of the two different cosmologies. We used the RLF 
model evaluated in Cosmology II of Jarvis et al. (2001c) to 
obtain the RLF used in Eq. 10. 

The best-fitting parameter logio -Do, the peak of the 
Gaussian distribution in log 10 of the projected linear size D 
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Figure 8. The estimated 151-MHz luminosity-redshift plane for 
the 6C** sample. The radio galaxies are represented by circles 
and the quasars by triangles. In the bottom-panel spectroscopic 
rcdshifts are taken into consideration (filled symbols). In both 
plots the star symbols represent the sources which were not de- 
tected in ii"-band (down to K ~ 21 mag in an 8-arcsec aperture), 
i.e. 6C**0737+5618 and 6C**0935+4348. The area between the 
horizontal lines is the region which contains the 'most luminous' 
sources according to our definition, and which is equivalent to the 
top-decade in luminosity defined by Jarvis et al. (2001c). These 
plots are for H M = 0.3, Q A = 0.7, H = 70 km s" 1 Mpc" 1 . 



(eq. 7 of Jarvis et al. 2001c), was also redefined due to the 
change in the set of cosmological parameters used. Since this 
is just a change in the Hubble constant, the relation between 
the linear sizes in the two cosmologies is then simply pro- 
portional to the ratio of their respective Hubble constants. 
For model C, logioDo becomes 1.97 in the cosmology used 
in this paper. All the other parameters remained the same 
(table 3 of Jarvis et al. 2001c). 



7.2.2 Using the best-fitting redshift estimates 

The redshift distribution of 6C** sources in the top-decade 
of luminosity, based on their best-fitting redshift estimates 
(.z os t) and on the spectroscopic redshifts available, is shown 
in Fig. 9. In general, the model redshift distribution is a 
fairly good approximation to the data. There is excellent 
agreement at z < 2, and it is also apparent that: (i) there is 
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Figure 9. The histogram shows the number of sources in the 
6C** sample with log 10 T151 > 27.61, drawn from the 'estimated 
+ spectroscopic' L151 — 2 diagram (Fig. 8: bottom-panel), binned 
in redshift with bin width Az = 1. The light shaded region repre- 
sents the sources which have spectroscopic redshifts. The solid line 
shows Jarvis et al. (2001c) model C predictions (rj = —0.06, con- 
stant co-moving space density) for the redshift distribution (Sec- 
tion 7.2.1). The dark shaded region represents 6C**0737+5618, 
for which the redshift estimate is a lower limit based on extrapo- 
lation of the K — 2 relation (see also the caption of Table 2). 

a slight deficit of sources predicted at 2 < z < 3; and, (ii) 
there is an excess of predicted sources at 2 > 3. However, we 
note that it is at these redshifts that the effects of the small 
number statistics become more important. Moreover, our 
redshift estimation method relies on the limited statistics 
of the 3CRR, 6CE, 7CRS and 6C* radio galaxies at z > 3, 
where the scatter in the K — z diagram is still poorly defined. 

Similar comments can be made when we use the RLF 
model to create an artificial L151 — z plane and compare 
it with the one estimated from the data, including the 
spectroscopic redshifts (Fig. 10). In general, it is apparent 
that the model predicts more sources (22) than the ones 
that are present in the estimated+spectroscopic diagram (17 
sources). However, we caution that, since both the estimated 
data and the simulations have independent Poisson errors, 
this comparison can only be regarded in a qualitative way. 

To explore the significance of these results, we apply the 
Kolmogorov-Smirnov test to compare the redshift distribu- 
tion derived from data and the predicted model distribution 
(Fig. 9). We find that they are not significantly different, 
with probability p — 0.29. Thus, we can conclude that the 
data are consistent with Model C of Jarvis et al. (2001c), 
the constant co-moving space density model. 

7.2.3 Taking the redshift probability distributions into 
account 

We now construct the redshift distribution of the top-decade 
sources by taking into account the redshift probability dis- 
tribution of each source (derived from the K — z diagram) , 
rather than just using the best-fitting redshift estimate. In 
practice, since the probability density functions (derived in 
Section 2.3) are normalised to unity, when the binned red- 
shift distribution is constructed each source contributes to 
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Figure 10. The top panel shows the estimated top-decade (de- 
limited by the dotted horizontal lines) of the L151 — 2 plane for 
the 6C** sample. The spectroscopic redshifts are included (by 
replacing estimated redshifts with spectroscopic ones when these 
are available), and are represented by the filled symbols. The star 
symbol represents 6C**0737+5618. The bottom panel shows a 
simulation of the top-decade of the L151 — 2 plane for a sample 
with the same flux limit, sky area and radio selection criteria as 
6C**. This simulation was generated using the best-fitting model 
C of Jarvis et al. (2001c), translated to our adopted cosmology. 



a number of redshift bins, with an amount proportional to 
the probability of the source lying within that bin. For each 
source, we calculate also the redshift which corresponds to 
a luminosity of log 10 L151 > 27.61, which is the lower limit 
of the top-decade, and only take into account contributions 
which lie above that redshift. All the sources are considered, 
i.e. we include also those sources whose best-fitting redshift 
estimate would put them below the lower-limit of the top- 
decade in luminosity, and consider the contribution from 
the high-redshift tail of their probability distributions. The 
spectroscopic redshifts, when available, are also taken into 
consideration. The sources with spectroscopic redshifts con- 
tribute to a single bin each, if their luminosity lies above the 
lower limit of the top-decade (see Table 3). 60**0737+5618, 
the source which is not detected in if-band and does not 
have a spectroscopic redshift, contributes also to a single 
bin (4 < z < 5). The resulting redshift distribution is pre- 
sented in Fig. 11. 

The distribution derived from the data appears to be 
in very good agreement with the model distribution. This 
is confirmed by the application of the chi-squared goodness- 
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Figure 11. The binned redshift distribution of sources in the top- 
decade (i.e. with log 10 L151 > 27.61), constructed by taking the 
redshift probability density functions of Section 2.3 into account, 
rather than just the best-fitting redshift estimates (see Fig 9), as 
described in Section 7.2.3. The bin width is Az = 1. Each source 
contributes to a number of bins, with an amount proportional to 
the probability of the source lying within that bin. The available 
spectroscopic redshifts are also taken into consideration. These 
contribute to a single bin each and are represented by the light 
shaded region. The solid line shows Jarvis et al. (2001c) model C 
predictions (n = —0.06, constant co-moving space density) for the 
redshift distribution (Section 7.2.1). The dashed, dotted and dot- 
dashed lines show the same model predictions but for declining 
co-moving space densities with n = —2.0, r) = —3.0 and r) = —4.0, 
respectively. The dark shaded region represents 6C**0737+5618, 
for which the redshift estimate is obtained based on extrapolation 
of the K — z relation, and a probability density function was not 
derived (see also the caption of Tabic 2). This source contributes 
to a single bin, the 4 < z < 5 bin. 



of-fit test, which gives x 2 = 1.2, with a probability p = 0.54 
of obtaining this, or greater than this value of \ 2 ■ 

We perform also the same analysis by omitting 
6C**0737+5618 in the 4 < z < 5 bin. This source is uniden- 
tified and it is debatable that its K-b&nd magnitude limit 
(K > 21 mag in a 8-arcsec aperture at the 3<r level) im- 
plies, by extrapolation of the K — z diagram, a very high 
redshift (z > 4). It is also possible that this source, like 
6C**0935+4348, is at a lower redshift (z ~ 2 - 3) but is 
perhaps at a very early stage of its formation and/or highly 
obscured, or is simply underluminous (see discussion in Sec- 
tion 4). However, a very high-redshift cannot be ruled out. 

Omitting 6C**0737+5618 from our analysis does not 
change the results significantly. We find \ 2 = 2.4, with p — 
0.36. Again, we conclude that the data are consistent with 
the constant co-moving space density model of Jarvis et al. 
(2001c). 

Finally, we compare the data with the predictions from 
five declining co-moving space density models, using Jarvis 
et al. (2001c) model C with: r\ — —2.0 (dashed line in 
Fig. 11), rj = -2.5, rj = -3.0 (dotted line in Fig. 11), 
q = —3.5 and rj = —4.0 (dot-dashed line in Fig. 11). The 
results of the chi-squared goodness-of-fit test are listed in 
Table 4. It can be seen that, although the data are also con- 
sistent with moderate declines by factors of 2 to 3, declines 



Omitting Notes 
6C**0737+5618 



Tj 


X 2 


P 


X 2 


V 




-0.06 


1.2 


0.54 


2.4 


0.36 


constant space density 


-2.0 


1.4 


0.24 


0.8 


0.36 


declining space density 


-2.5 


2.4 


0.12 


1.6 


0.21 


declining space density 


-3.0 


3.7 


0.05 


2.6 


0.11 


declining space density 


-3.5 


5.2 


0.02 


3.8 


0.05 


declining space density 


-4.0 


6.8 


0.009 


5.1 


0.02 


declining space density 



Table 4. The results of the chi-squared goodness-of-fit test to 
compare the model and data redshift distributions of Fig. 11. 
The first column lists the value of the parameter n, which is the 
power-law exponent describing the high-redshift co-moving space 
density (see Eq. 9) . Columns two and three give the values of chi- 
squarc x 2 , and the probability p of obtaining that value of \ 2 or 
greater, when including 6C**0737+5618 in the analysis. Columns 
four and five list the same values, when omitting this source. 

by a factor of 3.5 and 4.0 can be excluded at the ~ 2 — 3a 
level. 



7.2.4 The limitations of a filtered sample 

As discussed by Jarvis (2000), filtered samples such as 6C** 
have limitations when it comes to confirming a decline in co- 
moving space density. The major problem is that the lack 
of sources at a given redshift may not be due to a decline in 
their space density but to imperfections in the filtering tech- 
nique. For example, with respect to the 6C* sample, Jarvis 
(2000) estimated that the angular size selection is filtering 
out an increasingly large fraction of the sources with red- 
shift: from ~ 20% to ~ 30% between z — and z ~ 5; and 
~ 30 - 50% beyond z > 5. Jarvis (2000) concluded that 
samples with similar filtering criteria to that of 6C*, and in 
particular with flux-density limits similar to that of 6C**, 
are only able to confirm roughly constant or increasing co- 
moving space density at high redshifts. The presence of a 
decline would be difficult to interpret due to the uncertain- 
ties introduced by the filtering criteria. However, the filtering 
is helpful in placing strong lower limits on any decline. 



8 SUMMARY 

A method of redshift estimation, based on the K — z dia- 
gram of the 3CRR, 6CE, 6C* and 7CRS radio galaxies has 
been developed. Redshift probability density functions are 
derived for all of the 6C** sources which are identified with 
a near-infrared counterpart, i.e. for 66 of the 68 members of 
the sample. Comparison of the resulting redshift estimates 
with the subset of spectroscopic redshifts shows that our 
method is fairly robust whenever emission- line and/or non- 
stellar contributions to the /^-magnitudes can be neglected. 
The estimated redshift distribution has a median redshift of 
2 mcd ~ 1.6. However, we find that the quasars have their 
redshifts significantly under-estimated by our method. This 
is explained by the fact that the method is based on the K— z 
relation, which is only valid for radio galaxies. Removing the 
quasars from the distribution results in a median estimated 
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redshift of z « 1.7. This is similar to that of the 6C* sample 
(zmcd ~ 1.9) and is significantly higher than that of unfil- 
tered, complete surveys at the same flux density limit. We 
conclude that the filtering criteria were effective in biasing 
the 6C** sample to objects at high-redshift. 

The redshift distribution of the most luminous sources 
in the 6C** sample is compared with the predictions of the 
steep-spectrum RLF model of Jarvis et al. (2001c). We find 
that the 6C** data is consistent with a constant co-moving 
space density at z > 2.5, and moderate declines by factors of 
~ 4 can be excluded at the ~ 2 — 3a level. Although Jarvis 
ct al. (2001c) excluded these declines at the ~ 4a level, 
the additional data from 6C** provide an independent mea- 
sure. Thus, the two independent studies are in quantitative 
agreement with the result that any decline at high redshift 
is shallow. 

We note that our result is based on a redshift distri- 
bution which is uncertain for the following reasons: (i) a 
significant fraction of the sample is not identified spectro- 
scopically; (ii) the method of redshift estimation relies on the 
limited statistics of the 3CRR, 6CE, 7CRS and 6C* radio 
galaxies at z > 3 (the scatter in the K — z diagram at these 
high redshifts is still poorly defined); and (iii) the redshift 
estimates of quasars are systematically under-estimated by 
our method. Most of these are likely to lead to an under- 
estimate of the true median redshift of the complete sam- 
ple. Thus, although spectroscopically incomplete, with the 
6C** sample we have additional strong constraints on the 
high-redshift space density, with a sample that increases the 
number of powerful steep-spectrum sources, from complete 
samples at z > 2, by a factor of ~ 2. 

The work presented here could be significantly improved 
by obtaining spectroscopic redshifts for a larger fraction of 
sources in the 6C** sample. This would be particularly im- 
portant for the faintest sources (K > 19 mag), since these are 
the most probable z > 2 candidates. With spectra of these 
sources we should be able to obtain a tighter constrain on 
the co-moving space density at z > 2. 
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